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Gaussian quantum channels have recently attracted a growing interest, since they may 
lead to a tractable approach to the generally hard problem of evaluating quantum chan- 
nel capacities. However, the analysis performed so far has always been restricted to 
memoryless channels. Here, we consider the case of a bosonic Gaussian channel with 
memory, and show that the classical capacity can be significantly enhanced by employ- 
ing entangled input symbols instead of product symbols. 
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1. Introduction 

The evaluation of information capacities of quantum communication channels is of 
great interest in quantum information theory. In particular, an important question 
is to determine how much classical information can be processed asymptotically 
via a quantum channel. This problem has been solved, today, only for a few quan- 
tum channels, and it has been addressed only recently for bosonic channels, i.e., 
continuous-variable quantum channels acting on a bosonic field such as the elec- 
tromagnetic field ^. The classical capacity of a purely lossy bosonic channel was 
solved exactly very recently ^, while the case of noisy bosonic channels is already 
more involved. Actually, the classical capacity of the Gaussian bosonic channel, i.e., 
a continuous- variable quantum channel undergoing a Gaussian-distributed thermal 
noise, has been derived in although this result only holds provided that the opti- 
mal input ensemble is a tensor product of Gaussian states, as conjectured by sev- 
eral authors but not rigorously proven today (see e.g. ^ for recent progress on this 
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problem). All these studies, however, have been restricted to memoryless bosonic 
channels. 

The first analysis of the classical capacity of quantum channels with memory 
was performed in ^. In this work a qubit depolarizing channel with memory was 
investigated and the first evidence that entangled input states enhance the single- 
shot capacity was provided. It was shown that the single-shot capacity is always 
maximized cither by product states or by maximally entangled states at the input, 
and that the product versus entangled states optimal regime depends on a threshold 
value for the degree of channel memory. This threshold value turned out to be a 
function of the depolarization factor of the channel. The same features mentioned 
above for the depolarizing channel with memory were also later proved for a class 
of non- isotropic qubit channels with memory ^. 

In this paper, we investigate the capacity of a bosonic Gaussian channel that 
exhibits memory. We review the main results that were first presented in ^ and pro- 
vide more details about their derivation. We consider channels with a thermal noise 
that has a finite bandwidth. The resulting memory effect is modeled by assuming 
that the noise affecting two subsequent uses of the channel follows a bivariatc Gaus- 
sian distribution with a non- vanishing correlation coefficient, measuring the degree 
of memory of the channel. We prove that if the memory is non-zero and if the in- 
put energy is constrained, then the channel capacity can be significantly enhanced 
by using entangled symbols instead of product symbols. The relation between the 
degree of memory and the resulting optimal input entanglement is analyzed. 

The paper is organized as follows. In Sect. 2 we review the description and known 
results about memoryless bosonic Gaussian channels. In Sect. 3 we introduce the 
role of memory, and analyse the capacity by allowing the transmission of entangled 
input states. Finally, in Sect. 4 we summarise and discuss our results. 

2. Memoryless Bosonic Gaussian channels 

Let us define a memoryless bosonic Gaussian channel T acting on a mode of the 

electromagnetic field associated with the annihilation and creation operators a and 
a' , or, equivalently, the quadrature components q = {a + a^)/\/2 and p = i{a) — 
a)/\/2, satisfying the commutation relation [q,p] = i. If the input of the channel is 
initially in state p, we have 



where (fp = d^{/3) d'^{P) , while D(/3) = e'^"^-'^'" denotes the displacement op- 
erator (such that I a) = D{a)\0) with |0) being the vacuum state and |a) being a 
coherent state of mean value a). For a Gaussian channel, the kernel is a bivariate 
Gaussian distribution with variance N, namely 




(1) 



(2) 
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The channel then randomly displaces an input coherent state according to a Gaus- 
sian distribution, which results in a thermal state (N is the variance of the added 
noise on the quadrature components q and p, or, equivalently the number of ther- 
mal photons added by the channel). The Gaussian CP map effected by this channel 
can also be characterized via the covariance matrix. Restricting to Gaussian states 
with a vanishing mean value, a complete state characterization is provided by the 
covariance matrix 

The Gaussian channel can then be written as 

7-7+(o^j. (4) 

The coding theorem for quantum channels asserts that the one-shot classical 
capacity of a quantum channel T is given by 



Ci(T) = max 



S^p,T[p,]^-Y,p,S{T[p,]) 



(5) 



where S{p) = — Tr(plogp) is the von Neumann entropy of the density operator 
p. In Eq. (5), the maximum is taken over all probability distributions {pi} and 
collections of density operators {pi} satisfying the energy constraint 

^PiTr (piO^a) < n, (6) 

i 

with n being the maximum mean photon number at the input of the channel. For 
a monomodal bosonic Gaussian channel, it is conjectured that a Gaussian mixture 

of coherent states (i.e., a thermal state) achieves the channel capacity ^. The sum 
over i is replaced by an integral over a, where the input states 

p]^=\a}{a\ (7) 
are drawn from the probability density 

p{a) = ^e-^. (8) 
wn 

Thus, the one-shot classical capacity of the channel becomes 

Ci(T) =S{p)- j <fap{a) S [pT^) , (9) 
where we have defined the individual output states 



out 



and their mixture (saturating the energy constraint) 

1 

7r(n + N) 



d^ap(a)pr' = / d^f^e-^ |/3)(/3|. (11) 
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In order to calculate the entropy of a Gaussian state p, one computes the symplectic 
values of its covariance matrix 7, i.e., the solutions of the equation I7 — AJ| = 0, 
where 



i 

~i 



(12) 



It can be shown that these values always come as a pair ±A, so that the entropy is 
given by S{p) = g (|A| - i), where 

is the entropy of a thermal state with a mean photon number of x. Since the input 
states are coherent states with a covariance matrix 

y°^ifi?V (14) 



2 \0 1, 

the individual output states and their mixture p are associated with the co- 
variance matrices 

^ 2V 1 + 27V^' ^ 

l/l + 2(n + 7V) \ 
^ 2 V l + 2{n + N)J' ^ ' 

so that the one-shot capacity of the channel is 

C^iT) = g{n + N)-g{N). (17) 

3. Bosonic Gaussian channels with memory 
3.1. Multimodal Gaussian states 

Let us now generalize these notions to multimodal channels. Similarly to the 
monomodal case, s modes of the electromagnetic field, numbered from 1 to s, are 
associated to s pairs of annihilation and creation operators aj and at or, equiva- 
lently, of quadrature components qj and pj . Ordering these observables in a column 
vector 

rn 

R=[qi,pi,---,Qs,Ps] , (18) 

we define the mean vector m^*^ and the covariance matrix 7'*^ of an s-mode state 

p^**) as 

m^'*) = Tr p^'^lR, (19) 
7^^) =Tr{R- m^'^)p{R - m^''^)^ 'T^^Jh (20) 



I 



I 
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where each J, takes the form (12). We now focus on Gaussian states, which are 
completely specified by their mean vector and covariance matrix. The von Neuman 
entropy of such a state is then given by 

5(p^^^) = E -9(1^.1-^)' (21) 

where ±Xj arc the symplcctic eigenvalues of the covariance matrix 7^**) of the state; 
that is the solutions of the equation 

=0- (22) 

This will allow us to compute the one-shot classical capacity of bosonic Gaussian 
channels acting on s modes since this will still take the form (5), but with a 
generalized energy constraint 

i^EftTr(p(^)a]a,)<n. (23) 

j=l i 

3.2. Memoryless bimodal channel 

Consider two subsequent uses of a memoryless channel T, defining the bimodal 
channel 

T12M = J (ff3, (f(32 g(/3i,/32) 

X D(/3i) O D{(32) P D\(3i) O D\(32), (24) 

where 

(25) 

since the noise affecting the two uses is uncorrelated. Following Eqs. (18-20), we 
define the covariance matrix 712 of a bimodal state pu- We restrict ourselves to 
bimodal Gaussian states ^, characterized by 

7i2=(;r"), (26) 

where 71 is the covariance matrix associated with the reduced density operator 
pi — Tr2(pi2) of mode 1 (and similarly for 72), while (T12 characterizes the cor- 
relation and/or entanglement between the two modes. For a memoryless channel, 
the optimal input states are simply products of coherent states, with a covariance 
matrix 7J2 = 71^^ © 72^^ where 7}^ and 73^^ both take the form (14), while cr™ = 0. 
The optimal input modulation is a product of Gaussian distributions, 

, , 1 _|^il^+|»2£ , „, 

p{ai,a2) = ^^e « . 27 
It follows that the classical capacity of this channel is additive 5 x Ci(ri2) = Ci(T). 
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3.3. Bimodal channel with correlated noise 

Let us investigate what happens if the noise is correlated, for instance when the 
two uses are closely separated in time and the channel has a finite bandwidth. We 
assume that the noise distribution takes the general form 

,{0u02) = ^j=^e-^'^"-'', (28) 

where f3 = [3?(/3i), 3(/3i), 3?(/?2), 9(/?2)]'^ and 7^ is the covariance matrix of the 
noise quadratures, chosen to be 



7^ 



(29) 



N -xN \ 

N xN 
-xN N 
V ;riV N J 

Thus, the map T12 can be expressed by 712 >—> 712 + 7^, so that the noise terms 
added on the p quadratures of modes 1 and 2 are correlated Gaussians with vari- 
ance TV, while those added on the q quadratures arc anticorrelated Gaussians with 
variance N. In Sect. 3.4, we will see that the beneficial effect of entanglement dis- 
appears in a more symmetric noise model where the q and p noise quadratures are 
both correlated. The correlation coefficient x ranges from x = for a memoryless 
channel to a: = 1 for a channel with full memory. 

We now come to the central result of this paper. While we have seen that for a 
memoryless channel, the capacity is attained for product states, we will prove that 
for correlated thermal noise, the capacity is achieved if some appropriate degree of 
entanglement is injected at the input of the channel. Intuitively, if we take an EPR 
state, i.e., the common eigenstate of g+ = (gi + q2)/V^ and p_ = {pi — p2)/V^, 
it is clear that the noise on q+ and p- effected by the channel is reduced as x 
increases. This suggests that using entangled input states may decrease the effective 
noise, hence increase the capacity. However, EPR states have infinite energy so 
they violate the energy constraint. Instead, we may inject (finite-energy) two-mode 
vacuum squeezed states, whose covariance matrix is given by 

in in _ 1 / cosh2r 



7r^7?^i I. (30) 



in 1 / — sinh 2r 



'''' = 2[ sinh2rj' ^^^^ 

with r being the squeezing parameter. Note that purely classical correlations be- 
tween the quadratures in the input distribution p{ai,a2) also help increase the 
capacity when x > 0, so wc have to check that entanglement gives an extra en- 
hancement in addition to this. 

The mean photon number in each mode of the state characterized by Eqs. (30)- 
(31) is sinh^ r, so that the maximum allowed modulation (for a fixed maximum 
photon number n) decreases as entanglement increases. Remarkably, there is a 
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possible compromise between this reduction of modulation and the entanglement- 
induced noise reduction on g+ and p-. To show this, consider input states with 
sinh^ r = rjn; that is 

^1.2- 2 V l + 2r;n)' ^^^^ 



^in ^ 1 /-207n(r+^n) 
''^^ 2 V 2A/r;n(l + rjn) J ' 



(33) 



where r] measures the degree of entanglement and is used to interpolate between 
a product of vacuum states {rj = 0), which can be maximally modulated, and an 
entangled state {r] = 1), for which the entire energy is due to entanglement and 
no modulation can be applied. At the output of the channel, we get states with a 
covariance matrix 712^^ where 

out_ 1 /^l + 2?7n + 2^ \ 

- 2 ^ l + 277n + 27Vy' ' ^"^^^ 



^out 

^12 



1 f -2y/r]n{l+n) - 2xN 



2v^r;n(l +n) + 2xN 

while the mixture 712 of these states are characterized by 



(35) 



+[ (l-^)nj' 

assuming that the energy constraint is saturated. Here, y stands for the classical 
input correlation coefficient (to compensate for the noise, the q displacements need 
to be correlated, and the p displacements anti-correlated). 



3.4. Independent channels with phase-sensitive noise 

Before cvahiating the transmission rate achieved by these states, let us give some 
intuition about why entangled states could help enhance the transmission rate. 
More precisely, let us study this channel in a modified basis by considering the 
transformation 

(q+ = (qi + q2)/V2 
\p+ = ipi+P2)/V2 

which physically corresponds to a 50/50 beam splitter and maps the 1 and 2 modes 

to the + and — modes. Expressing the covariance matrix 7:^_ of the noise in this 
new basis confirms that the noise on the q+ and p- quadratures is reduced, while 



Q- = (?i - ?2)/v2 
P- = (Pi - P2)/V2 



(38) 
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the noise on the p+ and g_ is enhanced: 



/Nil 



7^ = 







V 



(39) 



\ 

N{1 + a;) 

N{l + x) 

N{l-x)J 

Moreover, in this basis, the channel introduces no correlation between the + and — 
modes, so that it acts as a product of two independent (memoryless) monomodal 
channels. These monomodal channels are phase-sensitive, in contrast to those acting 
on modes 1 and 2, since the variances of the q and p noise quadratures are different. 
As a consequence, their capacity is achieved by transmitting squeezed states ^, 
described by covariance matrices 

-2r 



7f = 



7- 











o2r 



(40) 



(41) 



where the squeezing parameter r is, as above, such that sinh r = rjn. Rotating 
back this product of squeezed states, corresponding to a diagonal covariance matrix 
7^ = 7^ © 7^, into the original basis, we obtain the two-mode squeezed states 
with covariance matrix 7}^ as defined in Eqs. (32)- (33). 

We may also continue studying this channel in the modified basis where all 
covariance matrices are diagonal. Indeed, the covariance matrix of states at the 
output of the channel is given by 7°li* = 7^1^ + 'y^_ . As for the mixture of output 
states, we have used rjn photons for squeezing, so that we are left with (1 — r])n 
photons for modulation. Distributing this photon number between the quadratures 
following the affecting noise then leads to a covariance matrix 



7+ 









ri)n 

{l-y){l-r])n 

(1 
















77)77. 

{l + y){l-n)nj 

where the parameter y allows to optimize the distribution of the modulation be- 
tween the quadratures. Up to the basis change, 7+_ is equivalent to the covariance 
matrix 712 characterized by Eqs. (34)-(37). 

Using this basis, we also easily understand why the transmission rate of a bosonic 
gaussian channel where the q and p noise quadratures are both correlated (or both 
anti-correlatcd) may not be enhanced by using entangled states. Indeed, the covari- 
ance matrix of the noise on the + and — modes takes in this case the form 



,(42) 



7t 



/N{l-x) 

Ar(i -x) 

iV(l-F 





\ 


x) 

N{l+x)J 



(43) 



I 



I 
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which is characteristic of a product of two independent phase-insensitive 
monomodal channels since the variance of the noise affecting both quadratures 

of a given mode is the same. The capacity of this channel is thus achieved by 
transmitting coherent states, which, after the change of basis, remain unentangled 
states. 

3.5. Entanglement- enhanced capacity 

In order to evaluate the transmission rate achieved by these states, we need first 
to compute the symplectic values A2*2^^ and A12 of 712^^ and 712, respectively (or, 

equivalently, of 7°^* and 7-1 ). The symplectic values ±Ai2 of a covariance matrix 

712 of the generic form (26) are the solutions of the equation 

|7i2-Ai2(JieJ2)|=0, (44) 

or, equivalently, the biquadratic equation 

At2-(l7i| + l72|+2|ai2|)AL + |7i2| =0. (45) 

Using Eqs. (34)- (37), we see that 7?2^* and 712 admit each one pair of doubly- 
degenerate symplectic values, namely 

^?2^* = ±V-out--out' A,2=±v/^^, (46) 

with 

«out = ^ + m + N, Wout = A/7?n(l + rin) + xN, 

u = \ +n + N, V = \/'qn{l + -qfi) + xN — y{l — rj)n. 

The transmission rate per mode is then given by 

i?(r;,y) =.9(|Ai2|-i)-.g(|A?]iV^) (47) 

We may now vary the degree of entanglement 77 and the classical correlation 
y to maximize the transmission rate R{r],y). In Fig. 1, where the level curves of 
R{ri,y) for n = 1 and = 1/3 are plotted, wc sec that when there is no memory 
{x = 0), the maximum lies at the origin {ri>,,y^:) = (0,0), while the transmission 
rate decreases as soon as we add classical correlation or entanglement. However, 
when a; > 0, the maximum is shifted towards non-zero r] and y values, that is, 
correlation and entanglement helps increase the transmission rate. For instance, 
when X = 0.7, the maximum is attained for (r;,, y,) = (0.19, 0.88). When the degree 
of memory x is further increased, the optimal transmission rate is achieved by 
using maximal classical correlation y = 1, e. g., for x = 0.9, we have {r]*,y*) = 
(0.30,1). Nonetheless, as x increases towards x = 1, the transmission rate may 
still be enhanced by increasing the degree of entanglement 77. At a; = 1, we have 
(r?*,y*) = (0.37,1). 

In Fig. 2, wc have plotted the transmission rate R optimized over the classical 
correlation y, as a function of the degree of entanglement rj. When x > 0, the 
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Fig. 1. Iso-transmission rate lines in the (r;, j/)-plane, where 77 is the degree of entanglement and 
y is the classical correlation coefBcient, for a thermal channel with a degree of memory a: = 0, 0.7, 
0.9 and 1. The mean number of photons is n = 1 at the input, while the number of added thermal 
photons is Af = 1/3. 



optimized rate i? over y increases with the degree of entanglement t] and attains a 
maximum at some optimal value r?* > 0. This confirms that entangled input states 
are useful to enhance the transmission rate provided that the channel exhibits some 
memory. 

It now suffices to maximize R with respect to both y and X] in order to find the 
channel capacity C (assuming that the conjecture * is verified and that no product 
but non-Gaussian states may outperform the Gaussian entangled states considered 
here). If we keep the signal-to- noise ratio njN constant, it is visible from Fig. 3 
that the optimal degree of entanglement 77* is the highest at some particular value 
of the mean input photon number n, and then decreases back to zero in the large-fi 
limit (except if x = or 1). Clearly, in this limit, the channel T tends to a couple of 
classical channels with Gaussian additive noise (one for each quadrature), so that 
entanglement cannot play a role any more . Fig. 4 shows the corresponding optimal 
value of the input correlation coefficient y* for the same values of the other relevant 
parameters. Note that, even in the classical limit n 00, some non-zero input 
correlation is useful to enhance the capacity of a Gaussian channel with a; > 0. 
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Fig. 2. Transmission rate R (maximized over the classical correlation coefficient j/) as a function 
of the input entanglement r] for a thermal channel with a degree of memory x. The mean number 
of photons is n = 1 at the input, while the number of added thermal photons is AT = 1/3. 
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Fig. .3. Optimal degree of input entanglement »7* as a function of the mean input photon number 
n for a fixed signal- to- noise ratio n/N = 3. 



4. Conclusion 

We have shown that entangled states can be used to enhance the classical capacity 
of a bosonic channel undergoing a thermal noise with memory. We determined 
the amount of entanglement that maximizes the information transmitted over the 
channel for a given input energy constraint (mean photon number per mode) and 
a given noise level (mean number of thermal photons per mode). For example, 
the capacity of a channel with a mean number of thermal photons of 1/3 and a 
correlation coefficient of 70% is enhanced by 10.8% if the mean photon number is 
1 and the two-modc> squeezing is 3.8 dB at the input. This capacity enhancement 
may seem paradoxical at first sight since using entangled signal states necessarily 
decreases the modulation variance for a fixed input energy, which seemingly lowers 
the capacity. However, due to the quantum correlations of entangled states, the noise 
affecting one mode can be partly compensated by the correlated noise affecting the 
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Fig. 4. Optimal degree of input correlation y* as a function of the mean input photon number 
n for a fixed signal-to-noise ratio n/N = 3. 



second mode, which globally reduces the effective noise. Interestingly, there exists a 

regime in which this latter effect dominates, resulting in a net enhancement of the 
amount of classical information transmitted per use of the channel. The capacity 
gain G, measuring the entanglement- induced capacity enhancement, is plotted in 
Fig. 5. It illustrates that a c;apacity enhancement of tens of percents is achievable 
by using entangled light beams with experimentally accessible levels of squeezing. 
It is interesting to notice that, unlike the qubit channels investigated so far, in the 
continuous variable case the capacity is always optimized by entangled states by 
varying the degree of entanglement, and that no threshold on the degree of memory 
is present. 
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